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FOREWORD 


This  work  was  performed  in  the  laboratories  of  the  Division  of 
Biological  Sciences,  National  Research  Council  of  Canada,  Ottawa,  Canada  as 
part  of  co-operative  DND-NRC  studies  of  the  biological  effects  of  nonionizing 
radiation. 


The  report  is  written  so  that  those  who  are  only  interested  in  using 
calibrated  probes  of  the  BRH  -  Narda  type  need  only  read  Chapter  1, 

Section  2.3  and  Chapters  5  and  6. 

The  remainder  of  the  report  will  be  of  interest  primarily  to  those 
who  are  involved  in  the  design,  construction,  evaluation  and  calibration  of 
implantable  miniature  electric-field  probes. 
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•— ~*A  set  i)f  tpiniature  probes  for  detennlning\electric  fields  In  tissue 
was  evaluated  and  calibrated  €^r  use  in  microwave  bibeffects  studies  at 
2450  MHz.  The  calibrations  in  &ir  and  tissue-equivalent  liquids  were  carried 
out  using  a  new  S-band  w^aypguide  technique.  The  air  calibration  using  wave¬ 
guide  has  an  accuracy  of  ^10%  compared  to^fcl8%  for  our\anechoic-chamber 
calibration.  The  average  probe  sensitivity  in  air  is  for  the 

five  probes  calibrated  and  varies  slightly  with  power  density  and  probe.  To 
estimate  probe  sensitivity  in  tissue,  a  section  of  waveguide  is  filled  with  a 
tissue-equivalent  liquid  and  is  separated  from  the  air-filled  waveguide  by  a 
very  thin  (0.25  mm)  planar  spacer.  The  probe  response  is  measured  as  a  function 
of  position  on  either  side  of  the  spacer  and  extrapolated  to  the  interface. 

The  ratio  of  probe  sensitivity  in  air  to  that  in  test  liquid  is  then  determined 
using  the  continuity  of  tangential  E  field  across  the  spacer.  Liquids  with 
dielectric  properties  simulating  both  wet  and  dry  tissues  were  qsed.  For  the 
water- glycerol  solution  modelling  wet  tissue  the  probes  are  3.0  4L0.6  times 
more  sensitive  thanfor  air.  When  used  in  tissue  the  total  calibration 

error  is  estimated  to  beX4>  30%  for'3^2^ancfs^0%  for  the  specific  absorption 
rate  at  the  site.  One  of  the  problems  with  this  probe  design,  the  mechanical 
weakness  near  the  %  ip ,  has  been  eliminated  in  newer  designs.  The  other  pro¬ 
blem,  that  the  size  of  the  probe  tip  (jk3  mm)  is  too  large  for  optimal  use  in 
small-animal  organs,  is  more  difficult  to  solve. 
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RESUME 


On  a  evalue  et  etalonnE  un  ensemble  de  sondes  miniatures,  utilises 
pour  determiner  les  champs  electriques  dans  des  tissus,  afin  d’Etudier  les 
effets  biologiques  des  micro-ondes  3  2450  MHz.  Lea  Etalonnages  dans  l’air  et 
dans  des  liquidea  Equivalent  a  des  tissus  ont  EtE  effectuEs  au  moyen  d’une 
nouvelle  technique  de  guide  d'ondes  de  la  bande  S.  La  prEcision  de  la 
calibration  dans  I'air  est  de  t  102  comparativeoent  a  ±  18Z  pour  l'etalonnage 
dans  notre  chambre  anechoique.  La  sensibilitE  moyenne  de  la  sonde  dans  l’air 
est  1,0  mV/mW  cm~2  pour  les  cinq  sondes  etudiees,  et  elle  varie  legerement  en 
fonction  de  la  densite  de  puissance  et  la  sonde  utilisEe.  Pour  ea timer  la 
sunsibllltE  de  la  sonde  dans  un  tissu,  on  remplit  une  section  de  guide  dfondes 
avec  un  liquide  Equivalent  3  un  tissu  et  on  Insole  du  guide  d’ondes  rempli 
d’air  avec  un  espaceur  plan  tres  mince  (0,25  mm).  La  rEponse  de  la  sonde  est 
mesuree  en  fonction  de  sa  position  de  cheque  cote  de  1 ’espaceur  et  extrapolEe 
a  lf interface.  Le  rapport  de  la  sensibilite  de  la  sonde  dans  l’air  a  celle 
dans  le  liquide  d’essai  est  ensuite  etabli  en  faisant  appel  3  la  continuitE 
du  champ  E  tangentiel  3  travers  1* espaceur.  On  a  employE  des  liquldes  ayant 
des  propriEtes  dielectriques  simulant  des  tissus  secs  et  hum! des •  Dans  la 
solution  £au-glycErol  reprEsentant  un  tissu  humide,  on  a  trouvE  que  les  sondes 
etalent  3,0  ±  0,6  fois  plus  sensiblea  3  E2  que  dans  l’air.  Employees  dans 
des  tissus,  on  a  estime  que  l’erreur  d’etalonnage  totale  des  sondes  est 
±  30Z  pour  E2  et  ±  402  pour  le  taux  d’absorption  specifique  a  1 ’emplacement  de 
mesure.  Un  des  deux  problemes  propres  a  cette  sonde,  soit  la  faible 
resistance  mecanique  pres  de  l’extremite,  a  etE  resolu  au  moyen  de  nouvel- 
les  formes.  L’autre  probletne,  c’est-a-dire  ext  remite  trop  grosse  (^3  mm) 
pour  etre  utilisEe  convenablement  dans  les  petits  organes  d’anlmaux,  est  plus 
difficile  a  resoudre. 
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1.  INTRODUCTION 

1 . 1  GENERAL 

For  microwave  bioeffacts  research  it  is  desirable  to  know  the  local 
rate  of  energy  deposition  at  the  site  of  action  of  the  microwave  radiation. 
The  rate  of  energy  deposition  is  usually  expressed  as  a  specific 
absorption  rate,  SAR  (W/kg) ,  and  is  known  from  the  basic  physics  to  be 
given  by 

SAR  -  (a/p)  E*  (1) 

where  Ej  is  the  total  electric  field-strength  (rms  V/m) ,  a^is  the 
electrical  conductivity  (mho/m)  and  P  is  the  density  (kg/m  )  of  the  medium. 
Since  o  and  p  are  both  known  quantities,  a  measure  of  E^  in  the  medium 
gives  the  SAR. 

Amongst  others,  Bassen  and  co-workers  (1,2)  have  produced  several 
versions  of  implantable  miniature  electric  field  sensors  and  have 
evaluated  them  for  the  cases  of  plane-wave  irradiation  in  air  and  in 
tissue-equivalent  models. 

1.2  GOALS  OF  THIS  STUDY 

We  have  obtained  seven  units  of  their  single-axis  dipole  design, 
produced  by  Narda  Microwave  Corporation.  The  aims  of  this  study  are: 

(i)  to  develop  and  evaluate  an  S -band-wave guide  calibration  technique  for 
a  2450-MHz  calibration  of  the  probes  for  use  both  in  air  and  in  tissue; 

(ii)  to  evaluate  the  performance  and  utility  of  this  particular  probe 
design,  and  (iii)  to  obtain  calibration  data  for  our  probes. 

1.3  PROBE  CONSTRUCTION 

The  probes  under  evaluation  were  produced  in  1978  by  Narda  Microwave 
Corporation  as  model  #25256.  The  design  was  that  of  the  (US)  Bureau  of 
Radiological  Health,  so  they  are  also  designated  as  BRH  model  10.  Our 
evaluation  was  of  five  original  probes,  numbers  5, 6, 7, 8  and  9,  and  of  two 
replacement  probes,  numbers  7B  and  8B. 

The  probe  design  is  shown  in  Fig.  1,  and  Fig.  2  gives  details  of 
the  dipole  antenna  construction  in  the  tip.  The  short  dipole  is 
constructed  of  two  lengths  of  gold  foil  with  a  total  length  of  1.3  mm  and 
has  a  low-barrier  Schottky  diode  as  the  detecting  element.  The  miniature 
antenna  is  encased  in  a  1.2-mm  thick,  2.3-mm  wide  planar  dielectric 
substrate  of  dielectric  constant  -  2  (2) .  The  rectified  dc  voltage  is 
taken  to  a  dc  voltmeter  through  the  two  lines  shown.  They  are  constructed 
of  high-resistance  material  to  prevent  direct  sensing  of  microwave  energy 
by  the  lines.  The  dipole  is  oriented  at  a  34°  angle  with  respect  to  the 
probe’s  long  axis  so  that  three  successive  rotations  of  the  probe  about 
its  own  axis  by  120°  results  in  dipole  orientations  in  three  mutually 
perpendicular  directions. 
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Fig.  2.  Design  details  for  the  miniature  dipole  antenna  structure  in  the 
probe  tip  (from  Bassenf  Franke,  Aphey  and  Aslan  (2),  with  permission). 
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The  probe  output  voltage,  V^,  for  probe  orientation  number  i 
(i  =  is  proportional  to  a  power  of  the  electric  field-strength, 

E^,  along  the  dipole  axis  for  that  orientation: 

"i  -  “k  <2> 

where  B^  is  an  empirical  constant  depending  on  the  medium  and  n  =  1  for 
an  ideal  diode,  0.95  ±  0.02  in  practice.  In  describing  our  results  we 
will  take  n  to  be  1  and  allow  Byj  to  vary  slightly  with  (or  V^)  . 

Summing  the  three  components  of  Eq.  (2)  gives  the  relationship 
between  the  total  voltage,  V-p,  and  the  total  electric  field.  Ex  : 

V  =  B  K2  (3) 

T  M  T 

1.4  PROBE  ELECTRICAL  AND  MECHANICAL  PROBLEMS 

At  first,  probe  voltage  readings  were  taken  with  a  high -input- 
impedance  (10^-0  ohms)  Keithley  Instruments  model  610C  solid-state 
electrometer.  An  RG  174/U  miniature  coaxial  cable  was  used  because  of 
its  light  weight  and  flexibility.  Cable  flexure  noise  was  ±0.2  mV, 
causing  no  problems. 

It  became  clear  over  a  period  of  six  months  that  the  use  of  the 
electrometer  was  occasionally  causing  the  miniature  probes  to  open  circuit, 
or  ’'burn-out"  during  the  process  of  switching  the  electrometer  off  one 
day  and  on  again  the  next.  An  investigation  of  electrometer  performance 
failed  to  find  any  evidence  of  significant  switching  transients  or 
abnormal  ac  or  dc  currents  from  the  electrometer  which  might  cause  this 
problem. 

H.  Bassen*  advised  us  that  one  of  his  BRH  model  10  probes  had  burned  out 
when  the  digital  voltmeter  it  was  connected  to  was  plugged  into  the  ac 
power  lines.  On  his  advice,  we  changed  to  using  a  Fluke  Model  8000A 
digital  multimeter  in  the  battery-operated  mode.  Probes  were  never  exposed 
to  microwave  radiation  when  not  connected  to  the  meter,  as  a  precaution 
against  possible  charge  accumulation  and  subsequent  discharge.  Under  these 
conditions,  no  more  probes  have  open-circuited.  Bassen*  believes  the  10-20“ 
M-ohm  input  impedance  of  the  voltmeter  is  what  eliminates  the  burn-out 
problem. 


In  addition  to  the  probe  burn-out  problem,  which  destroyed  three 
of  our  five  original  probes,  a  mechanical  weakness  was  found.  The  joint 
where  the  flat  2.5-cm-long  end  section  joins  the  3-mm-diameter  tube 
section  was  weak,  resulting  in  the  loss  of  one  probe  in  the  first  week  of 
use.  The  joint  was  strengthened  with  epoxy  glue  in  the  remaining  probes, 
which  were  used  with  as  much  care  as  possible,  but  despite  this,  one  more 
probe  broke  at  the  same  place.  Bassen  ert  al.  (2)  reported  that  several 
probe  tip  breakages  occurred  while  probes  were  pushed  through  a  sphere  of 
tissue-equivalent  material.  A  newer  version  of  the  probe,  incorporating 
an  epoxy-fiberglass  planar  tip,  has  eliminated  the  breakage  problem.* 

*  private  communication 
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Because  of  these  electrical  and  mechanical  problems  with  probes  as 
evaluation  tests  and  calibration  measurements  were  proceeding,  the  test 
results  are  incomplete  in  places.  In  particular,  no  useful  data  was 
obtained  from  probes  7  and  8  so  only  probes  5,6,7B,8B  and  9  appear  in  the 
tabulated  data.  At  the  time  of  writing,  only  (replacement)  probes  7B  and 
8B  are  still  operating. 

2.  CALIBRATION  IN  AIR-WAVEGUIDE  TECHNIQUE 

2.1  TECHNIQUE 

For  the  calibration  of  probe  response  in  air,  the  measurement 
system  of  Fig.  3  was  used  with  the  matched  load  as  the  final  waveguide 
component.  The  VSWR  of  either  test  cell  or  matched  load  was  measured 
with  a  commercial  slotted  line.  For  the  matched  load,  the  VSWR  is  less 
than  1^04.  The  probe  under  test  is  shown  in  its  holder  in  Figs.  1  and  3. 
The  54  angle  built  into  the  holder  allows  one  to  rotate  the  probe  to  the 
position  of  maximum  response  where  the  dipole  is  essentially  parallel 
(±  3  )  to  the  electric  field  in  the  waveguide.  Defining  the  voltage  for 
this  orientation  as  ,  itQwas  found  that  the  voltages  and  for  the 
other  two  orientations  120  from  the  first  were  small,  but  not  zero.  The 
average  result  for  all  probes  tested  was 

VT  =  V1  +  V2  +  V3  *  (1-07  1  °-02)  vr  (4) 

2 

The  calibration  of  vs  was  then  made  by  first  measuring  V^/V  for 
i  each  probe  at  one  power  level,  and  then  determining  V_  as  a  function  of 

T* 

The  microwave  frequency  was  determined  within  ±  3  MHz  with  a 
calibrated  coaxial  transmission  wavemeter.  The  frequency  varied  during 
operation,  and  from  day  to  day,  but  was  always  within  15  MHz  of  2450  MHz. 

Since  the  dipole  antenna  is  positioned  asymmetrically  inside  a 
dielectric  slab  it  is  possible  that  wave  distortion  effects  might  result 
in  a  difference  between  the  readings  for  a  wave  approaching  the  probe 
’’from  behind”,  as  in  Fig.  3,  as  compared  to  a  wave  approaching  from  the 
opposite  direction.  The  measured  difference  in  response  was  less  than 
2%  and  was  neglected. 

2.2  WAVEGUIDE  POWER  DENSITY  AND  ELECTRIC  FIELDS 

The  relationship  between  the  total  electric  field  at  the  center  of 
the  waveguide,  E^,  and  the  net  forward  power  in  the  waveguide,  Is  derived 
in  Appendix  I.  The  derivation  takes  into  account  the  fact  that  the  wave 
impedance  for  the  TE^q  mode  is  larger  than  377  ohms  by  the  ratio  of  the 
guide  wavelength  to  the  free-space  wavelength.  By  this  measurement  the 
value  of  BM  for  the  air  medium,  designated  BA,  was  determined.  Since  the 
probe  is  also  used  to  measure  the  power  density  in  free  space,  P^,  and 
since  P^  =  E^/377  (in  MRS  units)  the  equation  relating  VT  to  P^  may  also 
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be  written  as: 

VT  '  CA  V  <5> 

CA  depends  slightly  on  the  value  of  P^  (or  VT) .  Using  working  units  of 
mWem  ^  for  Pj,  mV  for  VT,  mV/mWcm  ^  for  CA  and  mV/V^m  2  for  b  it  can  be 
shown  that 

CA  =  3770  Ba  (6) 

2.3  RESULTS 

The  resulting  calibration  curves  of  v£  Pj  for  probes  7B  and  8B 
are  shown  as  Figs.  4  and  5,  respectively.  On  the  log-log  plot  shown,  the 
points  lie  on  a  straight  line  from  about  0.5  to  15  xrM/cwr  with  a  slope 
equal  to  the  power-law  index,  n.  The  probe  response  coefficient, 

CA  A  1  mV/mWcm  can  be  read  off  the  curve  at  different  values  of  Pi.  It 
has  been  tabulated  along  with  n  for  all  probes  tested  at  1,  10  and  50 
mW/cm^  in  Table  IV  of  section  3.3,  where  it  is  compared  to  other 
determinations  in  air. 

The  possible  errors  contributing  to  the  calibration  are  listed  and 
estimated  in  Table  I,  where  the  best  estimate  of  the  total  error  is  shown 
to  be  ±  10%. 

3.  CALIBRATION  IN  AIR  -  OTHER  TECHNIQUES 
3.1  ANECHOIC  CHAMBER  CALIBRATION 

Probes  7B  and  8B  were  also  calibrated  in  the  NRC  anechoic  chamber 
described  by  Assenheim  and  Hartsgrove  (3) .  It  was  found  impractical  to 
calibrate  the  probes  without  using  a  holder  for  setting  position  and 
orientation.  The  waveguide  holder  was  used,  and  to  average  out  the 
inevitable  field  perturbation  effects  of  the  holder,  three  different 
configurations  of  the  holder  with  respect  to  the  incident  plane  wave, 

|  defined  in  Table  II,  were  used. 

i 

The  probe,  in  its  holder,  was  placed  on  an  RF- transparent  block 

of  polystyrene  foam  at  a  distance  of  5.0  m  from  a  large  rectangular  horn 

antenna,  of  nominal  gain  22dB.  The  free-space  wavelength  for  2450  MHz  is 
A  -  12,2  cm  and  the  maximum  dimension  across  the  face  of  the  horn  is 
|  d  =  64  cm.  Using  these  figures,  the  distance  between  the  probe  and  the 

|  horn  was  calculated  to  be  1.49  d^/X,  which  puts  the  probe  essentially  in 

[  the  far-field  of  the  antenna. 

The  anechoic-chamber-callbration  results  are  presented  in  Table  III. 
Power-density  readings  performed  using  either  of  two  Narda  Microwave  model 
8305  radiation  monitors  or  a  Holaday  Industries  model  HI  1500  microwave 
survey  meter  agreed  within  5%.  The  absolute  accuracy  of  +  0.5  dB,  or 
!  12%,  stated  by  Narda  Microwave  in  their  manual  is  assumed  to  be  the 

main  error  in  power  density.  The  measured  probe  responses  in  Table  III 


* 


TABLE  I 
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Estimated  Errors  in  Waveguide 
Calibration  Procedure  at  2450  MHz 


Source  of  error 


Absolute  power  readings 
from  HP  435A  meter. 


Uncertainty  in  coupling 
coefficient  of  coaxial  reflectometer 
(±  0.2  dB) 


Variation  of  correction  factor, 
with  frequency, 

0.2% /MHz  over  ±  15  MHz, 


Vx> 


Reproducibility  of  the  ratio,  V^/V^, 

From  position  in  the  small 
standing  wave  pattern  (VSMR  <  1.05), 

Estimated  Total  Error 


Maximum  possible  (worst-case  sum) 

Root-mean-square  value, 

Best  estimate  (average  of 
the  preceeding  two  estimates) 


Estimated  error 
±~2% 

±  5% 

±  3% 

±  3% 

±  2% 

±15% 

±  7% 

±11% 


TABLE  H 

Definition  of  Configuration  used  for  Probe  and  Holder  in 
Anechoic  Chamber  Calibration  Measurements 


Configuration 

Probe  axis 

Probe  holder 
faces 

Probe  angle  for 
measurement  of 

(a) 

in  EK  plane^ 

horn 

dipole  parallel  to  E-field 

(b) 

in  EK  plane 

horn 

probe  rotated  60°  about  its 
own  axis  from  (a) 

(c) 

in  EH  plane 

sidewall 

as  for  (a) 

^Th e  e lectr i c  field,  E,  is  vertical.  The  magnetic  field,  H,  and  direction 
of  propagation  of  the  wave,  K,  are  horizontal. 


TABLE  III 


•4- 

CO 

CO 

o 

o 

O 

H 

+ 

A 

+i 

+i 

»—♦ 

NO 

rH 

> 

• 

• 

• 

a 

o 

rH 

00 

<r 

a* 

00 

on 

NO 

cd 

CO 

• 

• 

4-» 

>  o 

CO 

CM 

H 

CM 

O 

> 

\0 

NO 

On 

<u 

CM  • 

• 

• 

-Q 

>  CM 

N 

<r 

O 

r-i 

u 

CU 

NO 

ON 

vO 

i-H 

• 

• 

t>  vO 

o 

o 

co 

<y 

C  r-l 

o  jo 

•H  CTj 

<U  -u  H 
JO  cd 

o  u  c 

^  3  H 

pL,  00 

■H  t? 

4-1  <1) 

c  c 

O  1-1 
U  4-t 
OJ 
T3 


•CM 

^  a 

a)  u  u 

O  (0  S 
<34  C3  0 

<D  ^ 
T3 


rH  M 

oj  td  0) 
yO  *H  .O 

o  J-.  a 

M  <U  3 
Pm  W  3 


mean  ±  18%  :  ±  12%  from  P.  and  ±  6%  from  the  S.E. 

d 

mean  value  calculated  by  scaling  up  the  single  reading  in  configuration  (a)  in  proportion  to 
the  averages  determined  at  50  mW/cm^  for  the  three  configurations. 


12 


are  indicative  of  the  vertical  electric  field.  The  total  response,  VT,  is 
seen  to  depend  on  configuration.  The  best  estimate  of  V<p  is  taken  as  the 
mean  (±  S.E.)  of  the  three  configurations.  The  standard  error  of  6  or  7% 
in  VT  is  added  to  the  12%  power  density  error  to  give  an  estimated  total 
error  of  ±  18%  for  C^. 

3.2  NARDA  MICROWAVE  WAVEGUIDE  CALIBRATION 

Most  of  the  probes  were  supplied  from  Narda  Microwave  with  a  C^ 
value  determined  at  power  densities  of  1  and  10  raW/cm^.  Their  method 
involved  inserting  the  probe  through  a  small  hole  in  the  narrow  side  of  S- 
band  waveguide  and  positioning  it  so  that  the  dipole  sensed  the  full 
electric  field  at  the  center  of  the  waveguide.  Their  C^  values  were  given 
for  an  ideal  measuring  instrument  with  infinite  input  impedance  and  have 
been  reduced  12%  for  direct  comparison  to  our  data,  taken  with  a  10-Mohm- 
input-impedance  voltmeter. 

3.3  COMPARISON  OF  CALIBRATION  DATA  IN  AIR 

The  data  from  the  Narda  Microwave  waveguide  calibration,  our  own 
waveguide  determination  and  our  anechoic-chamber  measurements  are 
summarized  in  Table  IV.  The  C^  values  are  to  be  compared  keeping  in  mind 
the  estimated  uncertainties  of  ±  10%  for  our  waveguide  data  and  ±  18%  for 
our  anechoic-chamber  readings.  Narda  Microwave  did  not  estimate  their 
measurement  uncertainty.  The  data  are  generally  in  excellent  agreement. 
Our  values  for  the  power-law  index  lie  between  0.01  and  0.04  below  those 
of  Narda,  but  since  the  uncertainty  in  each  number  is  ±  0.02,  none  of  the 
differences  is  significant. 


WAVEGUIDE  CALIBRATION  IN  TISSUE- EQUIVALENT  LIQUIDS 


4.1  INTRODUCTION 

2 

The  probes  are  expected  to  produce  a  larger  voltage  for  the  same  E 
value  in  any  medium,  "M" ,  of  dielectric  constant,  c 1  >>  1,  than  for  air. 

Ln  other  words,  &m  =  BM  >  BA*  To  evaluate  &M  the  end  section  of 

waveguide  was  filled  with  diffeient  test  liquids  representing  a  range  of 
c*.  The  liquid-filled  waveguide  is  separated  from  the  air-filled  waveguide 
by  a  thin  planar  spacer  inserted  between  the  waveguide  flanges.  In  the 
air-filled  waveguide  only  the  TE^q  mode  exists  so  that  the  electric  field 
only  has  one  component,  which  is  tangential  to  the  surface  of  the  spacer, 
denoted  Et  (air).  For  our  probe  orientation,  V^  is  a  measure  of  Et  (air), 
i . e. ,  2 

W1  (air)  =  Ba  Et  (air) .  (7) 

Using  the  boundary  condition  that  the  tangential  component  of  E  is 
continuous  across  any  interface,  and  ignoring  (for  the  moment)  any  change 
in  E  within  the  thin  spacer,  the  tangential  component  of  E  in  liquid, 

K,  (liquid),  would  be  the  same  value  as  for  air.  This  should  be  true 
regardless  of  the  form  of  the  modes  propagating  in  the  liquid-filled 
waveguide.  That  is, 


Et  (liquid)  =  E  (air). 


(8) 


i 


Summary  of  Probe  Calibration  Data 
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The  response  in  the  liquid  medium  is  given  by 

„2 


V  (liquid) 


BM  (e')  Z"t  (liquid) 


(9) 


Using  Eq.  (8)  to  eliminate  the  E-fields  in  equations  (7)  and  (9)  we  get 


Bm  (c1)  V  (liquid,  at  interface) 
ba  V  (air,  at  interface) 


(10) 


By  measuring  as  a  function  of  position  on  each  side  of  the  spacer  and 
extrapolating  the  curves  up  to  the  interface,  the  ratio  of  probe  response 
coefficients  can  be  determined. 


The  electric  field  inside  the  test  liquid  is  known  (see  Appendix 
II)  to  decay  exponentially  due  to  the  energy  lost  in  heating  the  liquid. 
The  characteristic  decay  length  for  the  electric  field,  6  -  6  (e!  ,  o) 

is  defined  by  the  equation  1  1 

E(x)  -  E(o)  expt-x/fip  (11) 

where  E(o)  is  the  electric  field  at  the  interface  and  x  is  the  distance 
into  the  liquid  from  the  interface. 


Since  we  are  measuring  V  (x)  which  is  proportional  to  the  square 
of  E(x)  we  have  1 


V1(x)  =  V^(o) [exp(-x/61)  ]2 

=  V1(o)exp(-x/62)  (12) 

where  6^  =  6^/2  is  the  exponential  decay  length  for  (x) .  6^  has  also 

been  measured  for  the  different  liquids  for  comparison  to  the  value 
calculated  (in  Appendix  II)  from  the  known  electrical  properties  of  the 
liquids . 


4.2  TECHNIQUE 


An  attempt  was  made  to  use  a  thin-walled  tray  inside  the 
waveguide  to  contain  the  test  liquid.  The  results  were  found  to  depend 
critically  on  the  height  of  the  air  gap  between  the  top  of  the  liquid 
and  the  underside  of  the  waveguide.  As  a  result  of  this,  all  tests  were 
done  in  a  specially  constructed  cell  completely  filled  with  test  liquid. 
The  slot  in  the  top  of  the  cell  is  approximately  half  filled,  resulting 
in  a  power  density  radiating  out  of  the  slot  which  is  about  1%  of  the 
power  density  in  the  waveguide.  This  is  neglected  in  the  analysis. 

The  effect  of  spacer  thickness  on  the  measured  quantities 
%/bA  an<^  ^2  was  also  tested.  There  was  no  detectable  difference 
between  the  results  using  a  spacer  0.25  mm  thick  and  those  using  a 
1 . 37-mm-thick  spacer.  A  spacer  thickness  of  3.0  mm  altered  the  results. 
The  thinnest  spacer  (0.25  mm)  was  used  in  all  further  tests. 


An  analysis  of  preliminary  measurements  revealed  two  main  sources 
of  error  in  this  procedure  -  one  due  to  the  uncertainty  in  probe  position 
and  the  other  due  to  the  curve  extrapolation  process. 

To  eliminate  the  uncertainty  due  to  probe  position,  measurements 
were  taken  on  each  side  of  the  spacer  with  both  sides  air-filled.  Then 
the  position  of  each  curve  from  the  effective  center  of  the  spacer  was 
adjusted  for  the  best-fit  curve  through  the  data  on  both  sides  of  the 
interface.  This  procedure  reduces  the  position  uncertainty  to  ±  0.2  ram 
on  each  side  of  the  interface. 

The  other  appreciable  error  results  from  extrapolating  the  curves 
on  the  air  side.  The  combined  uncertainty  in  the  determination  of  B^/B^ 
resulting  from  these  two  errors  is  estimated  to  range  from  0.2  to  0.6. 
Approximately  2/3  of  the  error  comes  from  curve  extrapolation  which  in  turn 
depends  on  curve  shape  and  slope.  The  worst  case,  ±  0.6,  applies  to  the 
water-glycerol  solution  which  represents  wet  tissue.  A  computer  curve¬ 
fitting  program  would  help  here,  but  is  not  considered  to  be  worth  the 
effort . 

4 . 3  RESULTS 

Typical  results  for  V-^(x)  measured  on  both  sides  of  the  interface 
are  shown  for  probe  8B  in  Fig.  6.  As  described  above,  the  curves  (a)  for 
air  have  been  used  to  define  the  best  "zero"  for  distance  measurements 
from  the  center  of  the  spacer.  The  curves  (b) ,  (c)  and  (d)  using  the  same 
zero  point  show  the  increased  probe  response  in  test  liquids.  They  also 
show  clearly  the  exponential  decay  of  the  voltage  in  the  test  liquids, 
which  appears  as  a  good  straight  line  on  the  semi logarithmic  plot.  The 
curved  lines  on  the  air  side  of  the  interface  show  the  standing-wave 
pattern  in  air,  which  is  known  to  be  non-exponential. 

The  ratio  B^/B^  is  taken  from  the  curves  extrapolated  to  the  line 
representing  the  center  of  the  thin  spacer.  The  results  of  all  such 
measurements  performed  on  four  probes  using  four  test  liquids  are  presented 
in  Table  V.  Wet  tissue  can  be  seen  to  have  electrical  properties  closest 
to  the  glycerol-water  solution,  which  is  taken  to  be  the  best  model  of 
wet  tissue  in  this  study.  The  measurements  on  glycerol  and  water  were 
taken  to  help  understand  the  data  for  the  solution.  All  tests  were  done 
at  room  temperature,  19-22°C. 

All  measurements  of  6 2  are  presented  in  Table  VI.  The  data  show 
more  inconsistency  than  expected  from  the  estimated  uncertainty  of  about 
r  1  mm,  even  though  all  curves  in  liquid  were  good  exponentials.  Possible 
sources  of  variation  include  contamination  of  liquids  or  the  use  of  liquids 
not  at  room  temperature. 

4.4  DISCUSSION 

Regardless  of  the  ^  values,  the  B.,/B^  ratios  appc'ar  reliable 
within  the  error  of  ♦  0.2  to  0.6  stated  aoove.  The  average  response  ratio 
has  been  tabulated  for  each  test  liquid  and  each  probe.  There  are  probably 


-10  -5  0  5  10  15 

Distance  of  probe  dipole  from  centre  of  thin  spacer  (mm) 


Fig.  6,  Probe  voltage  vs  position  on  both  sides  of  a  thin  (0.25  mm)  space 
separating  air  from  four  different  test  fluids. 
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small  but  real  differences  between  some  of  the  liquids  and  probes,  as  seen 
in  the  table.  As  probes  7B  and  8B  were  replacement  probes  received  at 
different  times  after  the  original  set  of  probes,  their  fabrication  was 
slightly  different.  Both  had  dimensions  of  1.20  x  2.47  (i  0.03)  mm  for 
the  substrate  cross-section  compared  co  the  dimensions  1.20  x  2.35 
(i  0.03)  mm  for  all  of  the  original  set.  In  addition,  probe  7B  may  have 
been  encased  in  a  somewhat  different  substrate  from  the  others  since  it 
had  a  distinctively  greener  colour  than  the  others. 

The  test  solution  of  most  interest  is  the  water-glycerol  solution. 
The  average  response  ratio  is  3.00  ±  0.18  (SE)  .  Taking  into  account  the 
worst-case  bias  of  i  0.6  for  any  one  measurement  *  the  best  estimate  of  the 
error  in  the  ratio  is  *•  0.6.  The  response  of  probe  8B ,  used  for 
measurements  in  tissue,  can  best  be  estimated  by  multiplying  the  best  value 
for  Ujq/B^,  3.0  *  0.6,  by  the  ratio  of  the  average  response  of  probe  RB, 

2.15,  to  the  average  response  of  the  whole  group,  2.58,  i.e.,  the  best 

estimate  ol  B^/B^  for  probe  8B  for  the  glycerol-water  solution  is 
(3.0  1  0.6)  x  (2.15/2.58)  =  2.5  ±  0.5.  The  result  for  probe  7B  calculated 
i n  the  same  way ,  is  3.0  1  0.6. 

The  dependence  of  the  probe  response  ratio  on  dielectric  constant 
can  be  seen  in  Table  V.  The  probes  have  a  slightly  lower  response  in  the 
high-d i el ectr ic-cons tant  liquids  than  in  the  lower  two  ones.  This  is 
contrary  to  the  published  curves  of  Smith  (4),  which  were  calculated  for 
a  small  cylindrical  dipole  surrounded  by  a  cylinder  of  insulator.  In 
comparison,  we  have  a  small  planar  gold-foil  dipole  surrounded  by  a 
relatively  thin  planar  insulator.  Our  thin-insulator  results  are  more  in 
agreement  with  his  results  for  "sufficiently  thick11  insulation,  which  show 
that  for  all  e*  >  6  the  response  is  independent  of  ,  within  10%.  We 
measured  a  20  to  35%  (±  15%)  larger  response  in  liquids  of  lower  e*.  In 
view  of  the  moderate  increase  involved  and  the  difference  in  probe  geometry, 
the  disagreement  between  our  results  and  Smith’s  calculation  is  not 
considered  important.  The  most  important  conclusion  is  that  the  response 
ratio  is  so  nearly  independent  of  c1  that  the  probes  can  be  used  in  all 
types  of  wet  (and  some  dry)  tissue  without  knowing  e*  at  all. 

Bassen  et  al  .  (2)  recently  found  a  probe  response  ratio  of 

3.5  for  simulated  brain  material  and  5.4  for  simulated  muscle  material  in  a 
similar  probe  used  at  2.45  GHz.  They  also  reported  a  lower  response  ratio 
/or  materials  with  a  lower  t*  value.  As  the  probe  construction  was  not 
identical  and  the  dielectric  properties  have  not  yet  been  published,  no 
exact  comparison  can  he  made. 

Another  useful  conclusion,  deduced  from  the  curves  in  solution,  is 
that  there  is  no  alteration  in  probe  response  as  the  insulated  dipole 
approaches  as  close  as  1.7  mm  to  the  air-solution  interface,  as  predicted 
hv  Smith  (3). 


An  attempt  was  made*  to  determine  the  dominant  higher-order  mode 
in  each  liquid  by  comparing  the  measured  value  of  $2  with  the  values 
calculated  for  several.  TK[q  modes,  as  shown  in  Table  VI.  The  question 
marks  in  the  row  for  pure  glycerol  indicate  that  it  is  doubtful  iT  those 
modes  can  propagate  at  all.  Allowing  for  an  error  of  ±  10%  in  any  single 
measurement  of  ^2*  iL  cnn  he  seen  that  all  the  single  measurements  arc 
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62  (TKj^o  mode),  ns  expected. 

The  standard  error  in  6 2  for  water  is  0.9  mm.  From  Table  VII  in 
Appendix  II,  this  corresponds  to  a  temperature  change  of  2.2°C.  Each 
sample  of  liquid  was  used  at  room  temperature,  20.5  ±  1.5  °C,  but  sample 
temperature  was  not  measured  in  each  case.  It  is  clear  that  for  water, 
temperature  variation  could  account  for  much  of  the  scatter.  A  similar 
result  likely  applies  to  0.9%  saline  and  the  water-glycerol  solution. 

The  data  for  glycerol  may  indicate  a  TEiQmode  for  probes  7B  and 
8B  and  a  mixture  of  TE^q  and  TE30  modes  for  probes  5  and  6.  The  mode  of 
propagation  may  be  sensitive  to  conditions  at  the  interface.  For  example, 
if  the  thin  spacer  were  slightly  curved,  although  this  was  never  observed, 
it  might  excite  a  different  mode. 

Because  the  measured  standard  errors  are  comparable  to  the 
difference  in  62  between  successive  TE^jq  modes,  the  propagating  mode 
cannot  be  determined  for  certain.  It  is  most  likely  to  be  the  TE^q  mode 
since  that  mode  should  be  most  easily  excited  by  the  TE^q  mode  in  air  on 
the  other  side  of  the  thin  spacer. 

5 .  USE  OF  THE  PROBES  IN  DETERMINING  ELECTRIC  FIELD  AND  SPECIFIC  ABSORPTION 
RATE  IN  TISSUE 

5.1  RECOMMENDATIONS  FOR  THE  USE  OF  THE  PROBES 

Since  the  probes  break  easily  and  cannot  be  repaired  one  should 
take  care  in  all  aspects  of  handling  and  use.  As  the  miniature  cable 
connectors  are  awkward  to  use,  care  should  be  taken  when  connecting  the 
probes.  Caution  should  also  be  exercised  when  inserting  the  probes  Into 
tissue  and  in  rotating  them,  once  inserted.  They  should  not  be  used  in  an 
animal  that  could  possibly  move  and  break  them.  As  a  precaution  against 
damage,  it  is  a  good  idea  to  watch  the  probe  tip  during  all  handling 
processes . 

Electrical  burnout  problems  must  also  be  avoided  by  using  the 
probes  only  with  a  battery-operated  digital  voltmeter  with  an  input 
impedance  of  less  than  20  M  ohms.  As  a  precaution  against  charge 
accumulation  and  subsequent  discharge,  the  probe  tip  should  never  be 
exposed  to  microwave  radiation  when  it  is  not  connected  through  the  cable 
to  the  meter. 

5.2  DETERMINATION  OF  INTERNAL  ELECTRIC  FIELD  STRENGTH 

The  relationship  1  tween  total  voltage  and  total  electric  field, 
given  in  Eq.  3  above,  cai.  re-written  as 


(13) 
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Using  Eq.  6  to  eliminate  B^,  we  have 


2  (3770)  _^T 

T  ‘  (TO  '  CA 


(1M 


which  only  contains  the  measured  calibration  constants  CA  and 
If  CA  were  constant,  this  would  be  the  simplest  way  to  calculate  Ex* 
However,  since  CA  depends  slightly  on  Vp ,  we  can  obtain  the  most  useful 
formula  by  using  Eq.  5  to  eliminate  the  ratio  (VT/CA).  The  result  is 


E2  =  (  3770)  P  (V  ) 


<vv 
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where  Pj  (Vx)  is  the  power  density  determined  from  reading  the  air 
calibration  curve  for  the  probe  in  question  at  the  value  of  Vx  measured  in 
tissue . 


Probe  8B  was  used  in  all  tissue  measurements.  From  section  4.4, 
Bm/B/v  is  2.5  ±  0.5  for  this  probe,  giving  Ef  =  1508  Pd  *  Now  if 
Vx  =  10.0  mV,  for  example,  then  from  Fig.  5,  Pj  =  10.0  mW/cm^  and 
Ex  =  123  V/m.  The  error  in  is  i  20%  from  the  factor  of  2.5  ±  0,5  and 
±  10%  from  the  Pd(^x)  calibration  errors  (from  section  2.3)  for  a  total 
error  of  ±  30%.  This  assumes  perfect  contact  between  probe  and  tissue, 
or  the  physiological  saline  surrounding  the  tissue. 

5.3  DETERMINATION  OF  SPECIFIC  ABSORPTION  RATE  IN  TISSUE 


The  specific  absorption  rate  (SAR,  W/kg)  is  given  in  section 
1.1,  by  Eq.  1,  as 


SAR  =  (a/p) 

where  a  is  the  tissue  electrical  conductivity,  shown  in  Appendix  II  to  be  ^ 
1.9  i  0.2  mho/m  and  p  is  the  tissue  density,  taken  to  be  1  g/ml  =  10^  kg/m  , 
±  5%.  Combining  equations  1  and  15,  we  obtain 
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The  total  uncertainty  in  SAR,  including  i  30%  from  Ej,  i  10%  from  o  and 
t  5  7.  from  p,  is  about  i  *+0% . 


Continuing  the  previous  example  for  probe  8B  with  Bj.j/Ba  “  2,5  and 
Vx  =  1(1.0  mV,  the  SAR  la  29  +  11  mW/g,  including  the  total  error  of  *  40%, 
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5.4  MEASUREMENTS  OF  INTERNAL  ELECTRIC  FIELD  IN  MICE 

Cairnie  e^t  al .  (6)  have  reported  a  wide  range  of  dosimetry 
measurements  pertaining  to  the  exposure  of  mice  to  2450-MHz  radiation 
under  far-field  conditions  in  the  NRC  Anechoic  Chamber.  Part  of  that 
study  involved  internal  electric  field  measurements  in  the  testis  and 
abdomen  of  13  mice  in  5  different  orientations. 

One  major  problem  became  apparent  for  the  use  of  the  probes  in 
tissue  -  the  measurements  were  quite  inconsistent  between  mice  exposed 
under  the  same  conditions.  The  standard  deviations  of  the  measurements 
were  greater  than  half  the  mean  in  10/14  cases.  This  is  more  than  the 
usual  biological  variability. 

Two  explanations  for  the  large  scatter  in  the  results  are  possible. 
One  source  of  scatter  is  the  inconsistent  contact  between  probe  and  tissue. 
As  the  planar  probe  tip  is  rotated  in  one  testis  it  moves  the  tissue, 
changing  the  shape  of  the  testis  and  its  degree  of  contact  with  the  probe. 
Each  measurement  was  likely  done  with  a  different  area  of  probe  surface 
in  direct  contact  with  tissue,  as  opposed  to  small  air  pockets.  The  use 
of  physiological  saline  on  the  preparation  to  ensure  that  no  air 
contacted  the  probe  tip  had  limited  success,  partly  because  the  mouse 
could  not  usually  be  oriented  to  ensure  that  the  saline  did  not  flow  away. 

The  second  source  of  variation  might  result  from  the  very  strong 
dependence  of  readings  on  depth  of  penetration  into  the  tissue.  In  one 
brief  test,  where  the  reading  for  only  one  probe  angle  was  made,  the 
reading  of  =  7.8  mV  with  the  dipole  about  1  mm  into  testis  reduced  to 
a  constant  1,2  mV  at  3  or  4  mm  into  the  tissue.  This  2-  or  3-mm  decay 
distance  is  much  less  than  the  value  of  62  =  9.6  mm  calculated  for  the 
decay  of  E^  for  a  plane  wave  penetrating  a  flat  surface  of  material  with 
dielectric  properties  equal  to  testis.  The  effect  may  have  been  due  to 
the  curvature  of  the  surface  or  to  changes  in  probe  contact.  In  any  case, 
distance  into  the  testis,  which  may  be  crucial,  could  not  be  precisely 
controlled  in  those  experiments. 

6.  CONCLUSIONS  AND  RECOMMENDATIONS 

6.1  S-BAND  WAVEGUIDE  CALIBRATION  TECHNIQUE 

6.1.1.  Conclusions  for  Probe  Calibration  in  Air 

(i)  The  results  of  the  measurements  in  waveguide,  where  the  ratio 
E/H  ^  377  ohms,  must  be  converted  to  the  power  density  equivalent  to 
free  space  using  the  known  correction  factor,  A^/A.  This  correction 
is  accurate  and  easy  to  apply. 

(ii)  The  results  of  air  calibration  measurements  performed  in  three 
different  ways  by  two  groups  -  ourselves  and  Narda  Microwave  -  agreed 
within  the  experimental  errors. 

(iii)  Our  waveguide  technique  for  air  calibration  is  easier  than, 

and  more  accurate  than,  our  anechoic  chamber  technique  (±  10%  compared 
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to  ±  18%,  respectively).  The  difference  Is  not  important,  however,  since 
most  animal  experiments  involve  the  use  of  the  same  power  density  meters 
which  gave  rise  to  the  ±  18%  error  in  anechoic  chamber  calibrations  in 
the  first  place. 

6.1.2  Conclusions  for  Probe  Calibration  in  Tissue-Equivalent  Liquids 

(i)  The  probe  response  ratio,  B^/B^,  is  determined  by  a  direct  method 
involving  only  the  continuity  of  tangential  electric  field  at  an 
interface.  The  zero  of  position  on  each  side  of  the  interface  is  set 
by  measuring  a  standing  wave  pattern  with  air  on  both  sides  of  the 
interface . 

(ii)  Any  spacer  thickness  less  than  1.4  mm  is  likely  satisfactory. 

A  0.25-mm  spacer  was  used  in  all  tests. 

(iii)  Care  is  required  in  setting  and  measuring  position,  and  in  filling 
the  test  cell  to  the  correct  height  of  liquid. 

(iv)  The  method  is  very  well  suited  to  studying  response  as  a  function 
of  dielectric  constant. 

(v)  The  measurements  in  liquids  always  showed  an  exponential  decay, 
i.e.,  a  straight  line  on  semilog  paper.  The  measurements  in  air  were 
curves  on  semilog  paper,  giving  rise  to  appreciable  extrapolation  error. 

(vi)  The  measured  62  values  are  reasonable,  but  due  to  the  effects  of 
temperature  and  the  close  spacing  of  the  62  values  for  the  higher  order 
modes,  the  measured  6 2  cannot  be  used  to  precisely  determine  the  actual 
higher  order  mode  propagating.  The  TE^q  moc*e  *s  likely  the  only 
propagating  mode. 

6.1.3  Recommendations  for  Future  Work  in  Tissue-Equivalent  Liquids 

(.1)  Additional  accurate  data  on  the  dielectric  properties  of  pure 
liquids  and  solutions  is  required  at  2430  MHz. 

C i  1 )  Tests  in  liquids  at  higher  temperatures  c«^uld  be  used  to  check  the 
effect  of  tissue  temperatures  near  35  to  40  °C  on  probe  response. 

(iii)  For  further  work  to  identify  which  higher-order  modes  are 
propagating:  liquid  temperatures  must  be  controlled  and  measured;  a 

thicker  spacer  should  be  tried  to  ensure  no  bending  of  the  spacer 
interface;  and  the  VSWR  on  the  air  side  of  the  interface  should  always 
he  measured  as  a  check  on  the  consistency  of  the  propagating  mode  patterns, 

6.1.4  Comparison  with  the  Technique  of  Bassen  et  al» 

Bassen  et  al .  (2)  measured  fields  through  the  center  of  a  sphere 
of  tissue-equivalent  phantom  material  exposed  in  an  anechoic  chamber. 

The  following  comparisons  are  made: 


(i)  Both  techniques  require  the  construction  of  a  small  amount  of 
specilized  equipment. 
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(ii)  They  had  to  use  a  low-pass  0.5-Hz  filter  to  reduce  noise  resulting 
from  flexure  of  the  probe  as  it  was  moved  through  the  sphere.  We  had  no 
such  problem.  A  new  prototype  rigid  probe  design  has  significantly 
reduced  flexure  noise. 

(iii)  For  materials  equivalent  to  wet  tissue,  our  probe  response  ratio  of 
3.0  ±  0.6  agrees  with  their  value  of  3.6  for  brain-equivalent  material, 
but  disagrees  with  their  value  of  5.4  for  muscle-equivalent  material. 

(iv)  The  waveguide  and  sphere  techniques  are  complementary,  each  having 
its  own  advantages  and  disadvantages.  The  advantages  of  the  waveguide 
technique  are:  (1)  The  test  liquids  are  easier  to  work  with  than  the 
semi-solid  materials  used  in  the  sphere  method;  (2)  The  calculation  of 
the  field  pattern  on  both  sides  of  the  air-liquid  interface  is  much  simpler 
than  the  calculation  for  inside  the  sphere;  (3)  The  comparison  of  the 
measured  and  calculated  62  values  provides  a  check  on  the  accuracy  of  the 
measured  (or  assumed)  dielectric  properties;  (4)  The  errors  in  the 
method  can  easily  be  assessed. 

The  advantages  of  the  sphere  technique  are:  (1)  Spatial  averaging 
is  done  over  the  non-uniform  field  distribution;  (2)  The  extrapolation 
of  measured  curves  to  locations  inaccessible  to  the  probe  dipole  is 
avoided;  (3)  The  peak  E  field  is  located  by  passing  through  it. 

6.2  PERFORMANCE  EVALUATION  FOR  THE  BRH-NARDA  MINIATURE  ELECTRIC  FIELD  PROBE 
6.2.1  Performance  in  Air 


(i)  rrobe  output  voltage  has  been  found  to  be  proportional  to  E^n, 
from  0.5  to  15  mV,  with  a  slightly  decreasing  response  above  15  mV. 

The  value  of  n  is  0.95  ±  0.02,  very  close  to  unity. 

-2 

(ii)  Probe  sensitivity  ranges  from  0.89  to  1.31  mV/mWcm  for  the  five 
probes  tested.  The  difference  is  large  enough  that  a  separate  calibration 
curve  is  needed  for  each  probe. 

(iii)  The  small  probe  size  makes  it  useful  for  measuring  fields  in  any 
location  which  cannot  be  reached  by  conventional  power  density  meters 
because  of  their  large,  5-ctn -diameter,  sensors. 

6.2,2  Response  in  Tissue-Equivalent  Liquids 

(i)  The  increased  probe  response  in  tissue-equivalent  liquids  is 
expressed  as  the  "probe  response  ratio",  B^/B^.  This  was  found  to  be 
approximately  independent  of  the  liquid's  dielectric  constant,  within 
±  15%,  including  values  of  dielectric  constant  as  low  as  6.3  (for 
glycerol).  This  result  differs  from  Bassen's  (2)  and  also  does  not  agree 
with  Smith's  calculation  (5)  for  a  cylindrical  dipole  inside  an 
electrically  small  cylinder  of  insulator.  As  our  dipole  is  a  thin  foil 
embedded  in  a  thin  planar  insulator,  the  geometry  is  sufficiently 
different  that  another  result  might  be  expected.  Because  response  is 
independent  of  dielectric  constant,  the  probes  can  be  used  in  any  wet  or 
dry  tissue  where  they  can  be  inserted  without  breaking. 


! 
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(ii)  The  overall  average  value  of  B^/B^  measured  for  four  probes  in  four 
test  liquids  is  2.57  ±  0.11.  Including  possible  bias  in  the  zero  of 
position  measurements,  the  error  is  more  likely  about  +  0.4,  which  we 
consider  acceptable. 

(iii)  The  best  model  for  wet  tissue  is  the  water-glycerol  solution 
(70:30  v/v)  which  has  e*  =  58.5,  e11  -  19.5  and  B^/B^  =  3.00  f  0.18 
(mean  ±  SE  for  7  measurements  with  4  probes) . 

(iv)  The  response  ratio  averaged  over  all  liquids  was  found  to  be  almost 
independent  of  probe.  Only  probe  8B  (B^/B^  =  2.15  ±  0.14  (N  -  6)) 
differed  significantly  (t  =  2.698,  20  d.f.,  p  <  0.05)  in  its  response 
from  the  average  of  the  other  three  probes  (B^/B^  «  2.73  ±  0.12  (N  -  16)). 
Other  workers  using  their  own  set  of  similar  probes  need  not  calibrate 
them  in  tissue-equivalent  liquids  -  it  would  not  be  worth  the  effort. 

(v)  There  is  no  problem  using  the  probes  near  an  interface  between  any 
two  different  materials.  The  response  did  not  change  when  the  probe  was 
brought  as  close  as  1.7  mm  to  an  air-liquid  interface,  in  agreement  with 
Smith  ( 5) . 

6.2.3  Results  in  Tissue 

(i)  The  experience  of  Cairnie  ejt  al^.  (6)  making  measurements  in  mouse 
abdomen  and  testis  indicates  It  is  difficult  to  obtain  good  reproducibility 
of  measurements.  This  may  result  from  irregular  tissue  contact  or  a 
strong  dependence  of  results  on  probe  depth  in  tissue. 

(ii)  To  improve  the  consistency  of  tissue  contact,  the  mouse  should  always 
be  used  prone  so  that  applied  saline  will  not  run  off.  This  necessitates 
the  creation  of  an  electric  field  in  the  anechoic  chamber  of  either  vertical 
or  horizontal  polarization,  depending  on  the  mouse  orientation  being 

inves  tigated. 

(iii)  An  attempt  should  be  made  to  control  more  precisely  the  probe  depth 
and  position. 

6.3  RECOMMENDATIONS  FOR  IMPROVEMENTS  IN  PROBE  DESIGN 

(i)  A  more  rugged  mechanical  construction  for  the  probes  is  essential 
to  prevent  breakage.  This  has  been  done  for  the  newer  versions  of  the 
probe , 

(ii)  If  other  users  have  also  had  problems  with  electrical  burnout  of 
probes,  this  problem  must  be  investigated  and  solved.  It  does  not  seem 
likely  that  the  problem  lies  with  our  electrometer,  since  all  our  tests 
indicate  it  is  functioning  properly. 

(iii)  The  Microtech  EP-7S-1  six-pin  miniature  connector  has  been  found 
very  awkward  to  use.  There  is  too  great  a  danger  of  dropping  the  probe 
or  hitting  the  tip  against  something  while  struggling  to  connect  or 
disconnect  the  cable,  particularly  for  inexperienced  operators.  A  now 
2-pin  connector  should  be  selected. 
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(iv)  Rotating  the  flat  probe  tip  in  tissue  increases  the  chance  of  tip 
breakage.  It  also  may  change  the  degree  of  contact  between  probe  and 
tissue  which  would  greatly  affect  the  readings.  If  the  last  5  cm  of  the 
probe  were  made  into  a  rigid,  pointed  cylinder,  this  would  make  insertion 
easier  and  the  contact  more  consistent.  This  has  been  done  for  one  of 
the  newer  versions  of  the  probe. 


(v)  The  tip  diameter,  Just  less  than  3  mm,  is  too  large  for  use  in  many 
small-animal  organs,  such  as  the  testis  of  a  mouse. 
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APPENDIX  I 


CALCULATION  OF  POWER  DENSITY  IN  AIR-FILLED  WAVEGUIDE 


For  a  TE^q  mode  travelling  along  a  waveguide  with  a  (reflectionless) 
matched  load,  the  electric  (E)  and  magnetic  (H)  field  patterns  have  been 
described  by  Moreno  (7,  p  113).  The  (x,  y,  z)  co-ordinate  system  has  the 
origin  at  one  inside  corner  of  the  waveguide.  The  x-axis  runs  from  ,,0n  to 
"a"  across  the  broad  face  of  the  waveguide.  The  y-axis  runs  from  n0"  to  "b" 
across  the  narrow  face  of  the  waveguide,  and  the  z-axis  runs  the  length  of 
the  waveguide. 

The  E-  and  H-field  patterns  are  given  by 

£  =  E  sin(TTx/a)sin(u)t  -  Bz) 

y  o 

H  =  H  sin(7Tx/a)sin(u)t  -  Bz) 
x  o 

E  *  E  =H  =0 
x  z  y 

H  0  (but  is  not  needed)  . 


(17) 


Here,  oo  is  the  microwave  angular  frequency  (rad/sec),  B  is  the  propagation 
constant  and  EQ  and  Hq  are  constants  to  be  determined. 

The  net  power  flow  along  the  waveguide  may  be  calculated  from  the 
z  component  of  the  Poynting  vector; 


S  *  E  x  H 


(18) 


S  (x,y,t)  is  the  instantaneous  energy  flow  per  unit  area  (W/m  ). 

S  (x,y , t)  =EH  -EH  (19) 

z  x  y  y  x 

-  -E  H  ,  since  E  is  zero  everywhere.  The  minus 
sign,  which  just  indicates  the^d^rection  o£  energy  flow,  will  be  dropped, 
leaving 

S  (x,y,t>  =  E^H^  sin^(TTx/a)sin^(o)t  -  Bz).  (20) 

Defining  the  power  density,  P^(x,y)  as  the  time  average  of 
S  (x,y,t)  over  one  cycle  of  oscillation,  and  defining  P<j(c)  as  the  power 
density  at  the  center  of  the  waveguide,  the  above  equation  reduces  to 

Pd(x,y)  *  Pd  (c)sin2(7Tx/a)  ,  with  (21) 

PdU)  *  PdU  =  a/2,  y  =  b/2)  =  E  H  /2  (22) 

o  o 

=  E  H  (23) 

c  c 

where  Ec  and  Hc  are  the  rms  E-  and  H-field  strengths  at  the  center  of  the 
waveguide . 
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Defining  PQ  as  the  net  power  flow  along  the  waveguide,  it  is 
calculated  as  the  integral  of  the  power  density  over  the  cross-sectional 
area  of  the  waveguide  i.e., 

b  a 

P  “If  P ' (x,y)  dx  dy  (24) 

°  o  o  d 

=  a  b  (c)/2.  (25) 

Rearranging  equation  25  gives 


pi  (r)  =  _ jo  =  _o,  where  (26) 

d  V  '  ab  A 

A  =  ab  is  the  cross-sectional  area  of  the  waveguide. 

For  the  TF-^q  waveguide  mode  the  characteristic  wave  impedance  is 
not  the  same  as  Z  -  377  ohms,  which  applies  to  waves  in  free  space.  It  is 
given  by  Moreno  (/,  p  34)  as 

E  /H  =  Z  =  Z  (X  /A)  (27) 

y  x  o  g 

for  all  x  and  y.  X,  the  free-space  wavelength,  and  Ag,  the  guide  wavelength, 
are  related  to  each  other  and  the  cut-off  wavelength,  Ac,  by  the  well  known 
(Moreno  7,  p  123)  relationship 


where  X  =  2a. 
c 


(28) 


Applying  the  definition  of  Z 


to  the  power  density  formula  gives 


P'  (c)  =  E  H  =  E2/Z  =  ZH2. 
a  c  c  c  c 


2 

Since  the  probes  actually  measure  E  ,  not  P^,  and  since  Z  f  , 

the  calibration  in  waveguide  must  be  corrected  to  relate  it  to  the  power 

density  at  the  center  of  the  waveguide  equivalent  to  the  free-space 

situation,  F .  (c),  defined  using  Z  . 

a  o 


From  equations  27  and  39  we  have 


o 


Pd  (O 


=  P'  (O 
a 


(30) 
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The  formulas  will  now  be  evaluated  for  the  case  of  a  forward  power 
of  P  =  1  watt,  a  frequency  of  2450  MHz,  and  standard  S-band  waveguide 
(WR-284) . 

a  =  7.214  cm 

b  -  3.40  cm 

A  =  2a  -  14.43  cm 
c 

A  =  ab  =  24.51  cm^ 

P  =1  watt 
o 

A  =  12.24  cm 

A  =  23.09  cm 
g 

A  /A  =  1.89 
K 

fA  1 

Z  =  Z  =  (377  ohms)  (1‘89^  =  713  ohras 

o  (.A  j 

E  =  762  V/m  (rms) 
c 

H  =  1.070  A/m  (rms) 

c 

2  2 

?\  (c)  =  816  W/m  =  81.6  mW/cm 
d 

P  (c)  =  154  mW/cm*’ 
a 

Changing  the  frequency  by  10  MHz  to  2440  MHz  makes  Ag/A  =  1.91,  only 
17  different  from  1.89,  so  the  effect  of  small  changes  in  frequency  are 
ignored  in  this  work. 
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APPENDIX  II 


DIELECTRIC  PROPERTIES  AND  CALCULATED  DEPTH  OF  PENETRATION 
FOR  2450-MHz  RADIATION  IN  VARIOUS  MEDIA 


The  dielectric  properties  and  calculated  depth  of  penetration  t'oi 
several  different  modes  of  propagation  in  various  materials  of  interest  are 
listed  in  Table  VII.  The  goal  of  the  calculation  is  to  attempt  to  identify 
the  actual  mode  propagating  in  the  liquid  by  comparing  the  measured  .5.,  valu* 
to  the  calculated  ones  for  the  various  modes. 

The  standard  complex  number  notation  for  electromagnetic  waves 
propagating  in  any  homogeneous  isotropic  medium  is  used  in  this  Append i:-. 

The  complex  dielectric  constant,  or  relative  permittivity,  it,  therefore 
defined  as  C  =  c'  -  jt1'.  The  electrical  conductivity,  0,  is  given  by 
0  -  u)t  ?  ’  C0  y  where  eQ  =  b.85  picofarad/m,  u)  =  2ttI,  and  f  =  2.43  GHz. 

With  0  expressed  in  mho/m,  this  reduces  to  £ '  1  =  7.34  0. 

The  tabulated  dielectric  properties  are  taken  from  the  references 
given.  Since  and  0  are  both  in  common  usage,  both  are  given.  For 

physiological  saline  (9  g/1  NaCl),  the  data  of  Chou  and  Guy  ( M)  for  mammalian 
Ripger  solution  are  used,  as  the  two  solutions  are  electrically  very  similar. 
Testis  was  assumed  to  have  properties  equal  to  the  average  of  skeletal  muscle 
and  liver,  two  other  wet  tissues  that  have  been  measured  at  frequencies 
near  2.43  GHz.  In  particular,  the  muscle  and  liver  data  quoted  by  Hasted  (°, 
p  229)  for  1.00  and  3.00  GHz  were  used  in  a  linear  interpolation  calculation. 
The  dielectric  properties  of  water  were  calculated  using  linear  interpo:  t  ion 
in  temperature  and  frequency  from  the  data  of  three  different  papers  quo 1 1 
by  Hasted  (9,  pp.  4  3-43).  The  data  included  frequencies  of  1.744  ,  2.M  nr.d 
3.00  GHz  and  temperatures  of  20,  23,  30  and  40  °C. 

The  propagation  of  different  modes  in  waveguide,  when  described 
using  the  complex  number  notation,  has  E-  and  H- fields  which  vary  a  Lone  th. 
waveguide  as  exp(-^z).  y,  the  propagation  constant,  is  a  complex  number 
w  r  i  1 1  e  n  ,  1  s 


y  =  (i  +  j(i  t  •  1  > 

where  u  is  the  "attenuation  constant"  and  L  is  the  "phase  constant".  ;  and 
d  are  reaL  numbers  with  dimensions  of  inverse  length.  For  a  mode  propagating 
in  a  waveguide,  completely  filled  with  a  homogeneous,  isotropic  material  of 
complex  dielectric  constant  £ ,  ^  is  given  bv 

.  =  (2n/A)(0/\  )2-e)*  1  v- 

c 

where  \  is  the  cut-off  wavelength  for  the  mode  in  quest  ion  in  an  '»d 

waveguide.  Since  *  =  12.24  cm  this  reduced  to 

V  =  O.'.n  cm'1  ((12.24/X  )“  -  t. ’  +  Jl")*  (III 
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The  same  equation  applies  to  the  propagation  of  a  TEM  mode  in  an  unbounded 
medium  when  Ac  is  set  equal  to  infinity,  which  applies  to  the  actual 
situation  in  tissue. 
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Many  higher-order  modes  can  propagate  in  liquid-filled  waveguide. 
For  example,  for  £  =  49  -  jO,  /E1"  =  7  and  TE^  modes  as  high  as  the  TEqq 
mode  may  propagate.  In  our  case,  since  the  TE^q  ra0<^e  is  the  only  mode 
present  on  the  air  side  of  the  interface,  we  expect  most  of  the  energy  to 
go  into  the  TEmo  modes  in  solution,  but  not  necessarily  the  TE^q  mode. 

The  cut-off  wavelength  for  the  TEmo  mode  in  air-filled  waveguide,  denoted 
Acmo,  can  easily  be  shown  (Moreno  (7),  p  120) to  be  given  by 


A 

cmo 


2a  _  14.43  cm  # 
m  ro 


(34) 


Using  the  values  Aciq  =  14.43  cm,  Ac3q  =  4.81  cm  and  Ac5q  =  2.89  cm,  62  for 
each  of  the  corresponding  modes  in  each  test  liquid  was  calculated. 

63  for  the  TEM  mode  in  an  unbounded  medium  is  also  given  in  the 
table  to  show  that  it  never  differs  appreciably  from  the  value  for  the  TE^q 
mode . 

The  significant  effect  of  temperature  on  ($2  can  be  seen  in  the 
table  for  water,  where  the  rate  of  change  of  6 2  is  0.4  mm/deg  C. 
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